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We describe a protocol for continuously protecting un- 
known quantum states from decoherence that incorporates de- 
sign principles from both quantum error correction and quan- 
tum feedback control. Our protocol uses continuous measure- 
ments and Hamiltonian operations, which are weaker control 
tools than are typically assumed for quantum error correction. 
We develop a cost function appropriate for unknown quantum 
states and use it to optimize our state-estimate feedback. Us- 
ing Monte Carlo simulations, we study our protocol for the 
three-qubit bit-flip code in detail and demonstrate that it can 
improve the fidelity of quantum states beyond what is achiev- 
able using quantum error correction when the time between 
quantum error correction cycles is limited. 
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I. INTRODUCTION 

Long-lived coherent quantum states are essential for 
many quantum information science applications includ- 
ing quantum cryptography |^], quantum computation 
1^,^, and quantum teleportation [ij. Unfortunately, co- 
herent quantum states have extremely short lifetimes in 
realistic open quantum systems due to strong decoher- 
ing interactions with the environment. Overcoming this 
decoherence is the chief hurdle faced by experimenters 
studying quantum-limited systems. 

Quantum error correction is a "software solution" to 
this problem (^J^. It works by redundantly encod- 
ing quantum information across many quantum systems. 
The key to this approach is the use of measurements 
which reveal information about which errors have oc- 
curred and not about the encoded data. This feature 
is particularly useful for protecting the unknown quan- 
tum states that appear frequently in the course of quan- 
tum computations. The physical tools used in this ap- 
proach are projective von Neumman measurements that 
discretize errors onto a finite set and fast unitary gates 
that restore corrupted data. When combined with fault- 
tolerant techniques, and when all noise sources are below 



a critical value known as the accuracy threshold, quan- 
tum error correction enables quantum computations of 
arbitrary length with arbitrarily small output error, or 
so-called fault-tolerant quantum computation 0,^. 

Quantum feedback control is also sometimes used to 
combat decoherence [p[-pd|. This approach has the ad- 
vantage of working well even when control tools are 
limited. The information about the quantum state fed 
into the controller typically comes from continuous mea- 
surements and the operations the controller applies in 
response are typically bounded-strength Hamiltonians. 
The performance of the feedback may also be optimized 
relative to the resources that are available. For exam- 
ple, one can design a quantum feedback control scheme 
which minimizes the distance between a quantum state 
and its target subject to the constraint that all available 
controlling manipulations have bounded strengths jl^ . 

The availability of quantum error correction, which can 
protect unknown quantum states, and quantum feedback 
control, which uses weak measurements and slow con- 
trols, suggests that there might be a way to merge these 
approaches into a single technique with all of these fea- 
tures. Previous work to account for continuous time us- 
ing quantum error correction has focused on "automatic" 
recovery and has neglected the role of continuous mea- 
surement On the other hand, previous work on 
quantum state protection using quantum feedback con- 
trol has focused on protocols for known states and has 
not addressed the issue of protecting unknown quantum 
states however see for related work. 

The paper is organized as follows. In Sec. |l| we review 
quantum feedback control and introduce the formalism 
of stochastic master equations. In Sec. |l| we present the 
three-qubit bit-flip code as a simple example of a quan- 
tum error-correcting code and sketch the general theory 
using the stabilizer formalism. In Sec. [V we present our 



protocol for quantum error feedback control and derive 
an optimal non-Markovian feedback strategy for it. In 
Sec. we use Monte Carlo simulations to demonstrate 
this strategy's efficacy for the bit-flip code and compare 
it to discrete quantum error correction when the time be- 
tween quantum error correction cycles is finite. Section 
VI concludes. 
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II. QUANTUM FEEDBACK CONTROL 

A. Open quantum systems 

To describe quantum feedback control, we first need to 
describe uncontrolled open quantum system dynamics. 
Let S be an open quantum system weakly coupled to 
a reservoir TZ whose self-correlation time is much shorter 
than both the time scale of the system's dynamics and the 
time scale of the system-reservoir interaction. The Born- 
Markov approximation applies in this case and enables 
us to write down a master equation pO| describing the 
induced dynamics in S: 



p=-i [H, p] + ^I?[c^], 



(1) 



Here p denotes the reduced density matrix for S, H its 
Hamiltonian, and V a decohering Lindblad superopera- 
tor that takes a system-reservoir coupling operator (or 
jump operator) as an argument and acts on density ma- 
trices as 



I?[c]/3 — cpc^ 



If If 
—c'cp pc' c. 



(2) 



One way to derive this master equation is to imagine 
that the reservoir continuously measures the system but 
quickly forgets the outcomes because of rapid thermal- 
ization. The induced dynamics on S therefore appear 
as an average over all possible quantum trajectories that 
could have been recorded by the reservoir. 

What kinds of measurements can the reservoir continu- 
ously perform? The most general measurement quantum 
mechanics allows is a positive operator-valued measure 
(POVM) \Ej} acting on S. According to a theorem by 
Kraus 12^], the POVM {Ej} can always be decomposed 



as 



E 



(3) 



such that its stochastic action is p — > pj with probability 
Pj — tr {pEj ) , where 



Pj 



1 



tripE,) 



(4) 



This POVM is called a strong measurement when it can 
generate finite state changes and a weak measurement 
when it cannot POVMs that generate the master 

equation involve infinitesimal changes of state, and 
therefore are weak measurements. 

One reservoir POVM that results in the master equa- 
tion (|l|) is the continuous weak measurement with Kraus 
operators 



/ m 



.^c^ dt 



nfj,{dt) = VdtCf_,, 



1-1-= 



. , TO. 



(5) 
(6) 



Moreover, any POVM related to the one above via the 
unitary rotation 



(7) 



will generate the same master equation. We call each 
distinct POVM that generates the master equation when 
averaged over quantum trajectories an unravelling [ p3| of 
the master equation. 



B. Quantum feedback control 

The previous discussion of the master equation sug- 
gests a route for feedback control. If we replace the reser- 
voir with a device that records the measurement current, 
then we could feed the measurement record back into the 
system's dynamics by way of a controller. For example, 
the master equation (|I]) with m = 1 can be unravelled 
into the stochastic master equation (SME) ||20| , p^ 

dpc{t) = -i [H, p^{t)] dt 

+V[c]pcit)dt + H[c]pcit)dW{t) (8) 
dQ{t) = {c + c^)^dt + dWit), (9) 

where pc is the conditioned density matrix, conditioned 
on the outcomes of the measurement record Q{t), the 
expectation (a)^ means tr (pca), dW is a normally dis- 
tributed infinitesimal random variable with mean zero 
and variance dt (a Wiener increment p5[|), and 7i is a 
superoperator that takes a jump operator as an argument 
and acts on density matrices as 



H[c]p ^ cp + pc" — ptr \cp + pc^\. 



(10) 



This sort of unravelling occurs, for example, when one 
performs a continuous weak homodyne measurement of 
a field c by first mixing it with a classical local oscillator 
in a beamsplitter and then measuring the output beams 
with photodetectors |Q. The stochastic model (||-^) is 
flexible enough to incorporate other noise sources such as 
detector inefficiency, dark counts, time delays, and finite 
measurement bandwidth p^ ]. 

We can add feedback control by introducing a Q(t)- 
dependent Hamiltonian to the dynamics of pc- There are 
two well-studied ways of doing this. The first, and sim- 
plest, is to use Wiseman-Milburn feedback or cur- 
rent feedback, in which the feedback depends only on the 
instantaneous measurement current lQ{t) — dQ{t)/dt. 
For example, adding the Hamiltonian lQ{t)F to the SME 
(^) using current feedback leads to the dynamics 

dpc{t) = -i[H, pc{t)]dt 

+V[c]pc{t)dt + n[c]pc{t)dW{t) 
-i [F,cp^{t)+pcit)c^] dt 
+-D[F]p,{t)dt^t[F,p,{t)]dW 



dQ{t) = {c + c\dt + dW{t). 



(11) 
(12) 
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The second, and more general, way to add feedback is 
to modulate the Hamiltonian by a functional of the en- 
tire measurement record. An important class of this kind 
of feedback is estimate feedback [p^ , in which feedback 
is a function of the current conditioned state estimate 
Pc- This kind of feedback is of especial interest because 
of the quantum Bellman theorem [ p7| , which proves that 
the optimal feedback strategy will be a function only of 
conditioned state expectation values for a large class of 
physically reasonable cost functions. An example of such 
an estimate feedback control law analogous to the current 
feedback Hamiltonian used in (^l]) is to add the Hamil- 
tonian {lQ{t))cF — {c + c^)^F, which depends on what 
we expect the current Iq (t) should be given the previous 
measurement history rather than its actual instantaneous 
value. Adding this feedback to the SME (||) leads to the 
dynamics 

dpcit) ^-i[H, pc{t)]dt 

+V[c]pc{t)dt + H[c\pc{t)dW{t) 
^i{lQ),[F,Pc{t)]dt (13) 

dQ{t) = {c + c'')^dt + dW{t). (14) 



III. QUANTUM ERROR CORRECTION 

Although quantum feedback control has many merits, 
it has not been used to protect unknown quantum states 
from noise. Quantum error correction, however, is specif- 
ically designed to protect unknown quantum states; for 
this reason it has been an essential ingredient in the de- 
sign of quantum computers |p8|-^. The salient aspects 
of quantum error correction can already be seen in the 
three-qubit bit-flip code, even though it is not a fully 
quantum error correcting code. For that reason, we shall 
introduce quantum error correction with this example 
and discuss its generalization using the stabilizer formal- 
ism. 



A. The bit-flip code 

The bit-flip code protects a single two-state quantum 
system, or qubit, from bit-flipping errors by mapping it 
onto the state of three qubits: 

|0) ^ 1000) ^ |0) (15) 
|1)^|111)^|I). (16) 

The states |0) and |i) are called the basis states for 
the code and the space spanned by them is called the 
codespace, whose elements are called codewords. 

After the qubits are subjected to noise, quantum error 
correction proceeds in two steps. First, the parities of 



neighboring qubits are projectively measured. These are 
the observables[] 

Mo = ZZI (17) 
Ml = IZZ. (18) 

The error syndrome is the pair of eigenvalues (mo, mi) 
returned by this measurement. 

Once the error syndrome is known, the second step is to 
apply one of the following unitary operations conditioned 
on the error syndrome: 



(--1,+1)^A// (19) 

(-1,-1) ^/X/ (20) 

(+1,-1) ^//X (21) 

(+1,+1)^///. (22) 



This procedure has two particularly appealing char- 
acteristics: the error syndrome measurement does not 
distinguish between the codewords, and the projective 
nature of the measurement discretizes all possible quan- 
tum errors onto a finite set. These properties hold for 
general stabilizer codes as well. 

If the bit-flipping errors arise from reservoir-induced 
decoherence, then prior to quantum error correction the 
qubits evolve via the master equation 

dpnoisc = l{V[XII] + V[IXI] + V[IIX])p dt, (23) 

where ^dt is the probability of a bit-flip error on each 
qubit per time interval -(- dt]. This master equation 
has the solution 

a (t) po 

+b (t) (XllpaXII + IXIpoIXI + II X poll X) 
+c (t) [XXIpoXXI + XIXpoXIX + IXXpoIIX) 



+d{t)XXXpoXXX, (24) 
where 

a{t) = (1 + Se-^-*-* + 3e-4T* + e-^T*) /8 (25) 

b{t) = (1 + e-27t _ g-^it _ g-67t-) /g (26) 

c(<) = (1 - e-'^* - e-^''* + e-^T*) /8 (27) 

d{t) = (l - 3e-2^* + 3e-^T* - 6"^^*) /8. (28) 



The functions a{t)~d(t) express the probability that 
the system is left in a state that can be reached by zero, 
one, two, or three bit-flips from the initial state, respec- 
tively. After quantum error correction is performed, sin- 
gle errors are identified correctly but double and triple 



^ We use the notation of |2^] in which X, Y, and Z denote the 
Pauli matrices ax, cjy and respectively, and concatenation 
denotes a tensor product (e.g., ZZI = (gi I). 



3 



errors are not. As a result, the recovered state, averaged 
over all possible measurement syndromes, is 

p={a{t)+b (t)) po + (c (t) + d (t)) XXXpoXXX. (29) 

The overlap of this state with the initial state depends 
on the initial state, but is at least as large as when the 
initial state is |0); namely, it is at least as large as 

Fg = (2 + 3e^27t _ g-67t^ /4 

~ 1-3(70'- (30) 

Recalling that a single qubit subject to this decoherence 
has error probability p — ^t, we see that, when applied 
sufhciently often, the bit-flip code reduces the error prob- 
abihty on each qubit from 0{p) to 0{p^). 

B. Stabilizer formalism 

The bit-flip code is one of many quantum error cor- 
recting codes that can be described by the stabilizer for- 
malism [Q. Let C be a 2'^-dimensional subspace of a 
2"-dimensional n-qubit Hilbert space. Then the system 
can be thought of as encoding k qubits in n, where the 
codewords are elements of C. Let us further define the 
Pauli group to be P„ = {±l,±i} (E) and 
let the weight of an operator in Pn be the number of 
non-identity components it has when written as a tensor 
product of operators in Pi. The stabilizer of C, S{C), 
is the group of operators which fix all codewords in C. 
We call C an [[n, fc, d]] stabilizer code when (a) the n — k 
generators of S{C) form a subgroup of P„ and (6) d is the 
smallest weight of an element in P„ \ S{C) that commutes 
with every element in S{C). 

In this general setting, quantum error correction pro- 
ceeds in two steps. First, one projectively measures the 
stabilizer generators to infer the error syndrome. Sec- 
ond, one applies a unitary recovery operator conditioned 
on the error syndrome. The strong measurement used 
in this procedure guarantees that all errors, even unitary 
errors, are discretized onto a finite set. For this reason we 
will sometimes refer to this procedure as discrete quan- 
tum error correction. When the noise rate is low and 
when correction is applied sufficiently often, this proce- 
dure reduces the error probability from 0{p) to 0{p^). 

IV. CONTINUOUS QUANTUM ERROR 
CORRECTION VIA QUANTUM FEEDBACK 
CONTROL 

In this section, we present a method for continuously 
protecting an unknown quantum state using weak mea- 
surement, state estimation, and Hamiltonian correction. 
As in the previous section, we introduce this method via 
the bit-flip code and then generalize. 



A. Bit-flip code: Theoretical model 

Suppose p is subjected to bit-flipping decoherence as 
in (|23|); to protect against such decoherence, we have 
seen that we can encode p using the bit-flip code 
|r^ . Here we shall define a similar protocol that operates 
continuously and uses only weak measurements and slow 
corrections. 

The first part of our protocol is to weakly measure the 
stabilizer generators ZZI and IZZ for the bit-flip code, 
even though these measurements will not completely col- 
lapse the errors. To localize the errors even further, we 
also measure the remaining nontrivial stabilizer operator 
The second part of our protocol is to apply the 
slow Hamiltonian corrections XII, IXI, and IIX corre- 
sponding to the unitary corrections XII, IXI, and IIX, 
with control parameters Afe that are to be determined. If 
we parameterize the measurement strength by k and per- 
form the measurements using the unravelling (H-^, the 
SME describing our protocol is 

dpc = -f{V[XII] + V[IXI] + V[IIX])pcdt 
+k{V[ZZI] + V[IZZ] + V[ZIZ])pcdt 

+^{n[zzi]dWi + n[izz]dW2 
+ n[ziz]dW3)pc 

~i[F,p,]dt (31) 
dQi = 2k{ZZI)^ dt + ^dWi (32) 
dQ2 = 2k{IZZ\ dt + y/^dW2 (33) 
dQs = 2K{ZIZ\dt + (34) 

where 

F ^ XiXII + X2IXI + X3IIX (35) 
is the feedback Hamiltonian having control parameters 
. . 

Following the logic of quantum error correction, it is 
natural to choose the Afe to be functions of the error syn- 
drome. For example, the choice 

Ai = ^(1 - {ZZI)J{1 + (/ZZ)J(1 - 

A2 = ^(1 - {ZZI),){1 - {IZZI){1 + {ZIZD 

A3 = ^(1 + (ZZ/)J(1 - {IZZ),){1 - {ZIZ)J, (36) 

where A is the maximum feedback strength that can be 
applied, is reasonable^: it acts trivially when the state is 



^The modest improvement gained by this extra measure- 
ment is offset by an unfavorable scaling in the number 
of extra measurements required when applied to general 
[[n, k, d\] codes having 2""*^ stabilizer elements and only n — k 
generators. 

■^The factor of | is included to limit the maximal strength 
of any parameter to A. 
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in the codespace and applies a maximal correction when 
the state is orthogonal to the codespace. Unfortunately 
this feedback is sometimes harmful when it need not be. 
For example, when the controller receives no measure- 
ment inputs (i.e., k = 0), it still adds an extra coherent 
evolution which, on average, will drive the state of the 
system away from the state we wish to protect. 

This weakness of the feedback strategy suggests that 
we should choose our feedback more carefully. To do 
this, we introduce a cost function describing how far away 
our state is from its target and choose a control which 
minimizes this cost. The difficulty is that our target is 
an unknown quantum state. However, we can choose 
the target to be the codespace, which we do know. We 
choose our cost function, therefore, to be the norm of the 
component of the state outside the codespace. Since the 
codespace projector is Uq = \{III + ZZI + ZIZ + IZZ), 
the cost function is 1 — /, where f{p) — tr{pllc)- Under 
the SME ( |3l| ) , the time evolution of / due to the feedback 
Hamiltonian F is 

2Xi{YZI + YIZ)^ 
+2\2{ZYI + IYZ)^ 
+2X3{ZIY + IZYI. (37) 

Maximizing ffb minimizes the cost, yielding the optimal 
feedback coefficients 

Ai = Xsgn{YZI + YIZ)^ 

\2 = Asgn(zr/ + /rz), 

A3 = A sgn(Z/r + IZY)^ , (38) 

where, again, A is the maximum feedback strength that 
can be applied. 

This feedback scheme is a hang-hang control scheme, 
meaning that the control parameters Afc are always at 
the maximum or minimum value possible (A or —A, re- 
spectively), which is a typical control solution both classi- 
cally ||3l[] and quantum mechanically In practice, the 
bang-bang optimal controls ( ^8| ) can be approximated by 
a bandwidth-limited sigmoid, such as a hyperbolic tan- 
gent function. 

The control solution (^) requires the controller to in- 
tegrate the SME (|l]) using the measurement currents 
Qi{t) and the initial condition pc- However, typically the 
initial state Pc(0) will be unknown. Fortunately the cal- 
culation of the feedback (|3|) does not depend on where 
the initial condition is within the codespace, so the con- 
troller may assume the maximally mixed initial condition 
Pe = 5(|5)(0| + |1)(1|) for its calculations. This property 
generalizes for a wide class of stabilizer codes, as we prove 
in the appendix, and we conjecture that this property 
holds for all stabilizer codes. 



B. Intuitive one-qubit picture 

Before generalizing our procedure, it is helpful to gain 
some intuition about how it works by considering an 
even simpler "code": the spin-up state (i.e., |0)) of a 
single qubit. The stabilizer is Mq = Z^ the noise it pro- 
tects against is bit flips X, and the correction Hamil- 
tonian is proportional to X. The optimal feedback, 
by a similar analysis to that for the bit-flip code, is 
F — Asgn(y)cX, and the resulting stochastic master 
equation can be rewritten as a set of Bloch sphere equa- 
tions as follows: 

d{X)^ ^ -2k{X)^ dt - 2 (Z)^ dW (39) 
d{Yl - -2j{Y),dt ~ 2k{Y), - 2^{Y), {Z),dW 

-2X{sgn{Yl){Z),dt (40) 
d{Z)^ = -2j{Z)^dt + 2V^(1 - {Z)^)dW 

+2A(sgn(r)J(r),dt. (41) 

The Bloch vector representation ((X), (F), (Z)) [|| of 
the qubit provides a simple geometric picture of how it 
evolves. Decoherence (the 7 term) shrinks the Bloch vec- 
tor, measurement (the n terms) lengthens the Bloch vec- 
tor and moves it closer to the 2;-axis, and correction (the 
A term) rotates the Bloch vector in the y-z plane. Fig. |l| 
depicts this evolution: depending on whether the Bloch 
vector is in the hemisphere with (Y) > or (Y) < 0, 
the feedback will rotate the vector as quickly as possi- 
ble in such a way that it is always moving towards the 
codespace (spin- up state). Note that if the Bloch vector 
lies exactly on the z-axis with (Z) < 0, rotating it either 
way will move it towards the spin-up state — the two di- 
rections are equivalent, and it suffices to choose one of 
them arbitrarily. 



z=l 




z=-l 
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FIG. 1. Bloch sphere showing the action of our feedback 
scheme on one qubit. Wherever the Bloch vector is in the 
y-z plane, the feedback forces it back to the spin-up state, 
which is the codespace of this system. All the vectors shown 
lie, without loss of generality, in the a; = plane. 



C. Feedback for a general code 

Our approach generalizes for a full [[71, fc, d]] quantum 
error-correcting code, which can protect against depolar- 
izing noise |^ acting on each qubit independently. This 
noise channel, unlike the bit-flip channel, generates a full 
range of quantum errors — it applies either X, Y, or Z 
to each qubit equiprobably at a rate 7. We weakly mea- 
sure the n — k stabilizer generators {Mi} with strength 
K. For each syndrome m, we apply a slow Hamiltonian 
correction Fm with control strength Am, the weight of 
each correction being d or less. The SME describing this 
process is 

n n — k 

(^Pc= 7 ^ ^{V[af'>])pcdt + Kj2'^[Mi]Pcdt 

j—x.y,z i—1 l—l 
n — k R 

7^[Mi]dl¥,p, -iY,MFr,Pc]dt. (42) 

The number of feedback terms R needed will be less 
than or equal to the number of errors the code corrects 
against. The reason that this equality is not strict is that 
quantum error correcting codes can be degenerate, mean- 
ing that there can exist inequivalent errors that have the 
same effect on the state — a purely quantum mechanical 



property |28|. 

We optimize the Xr relative to a cost function equal 
to the state's overlap with the codespace. For a general 
stabilizer code C, the codespace projector is 



n, 



1 



c — 



2"- 



1=1 



and the rate of change of the codespace overlap due to 
feedback is 



ffb = -ltrJ2Xr[Ilc,Fr]f 



r=0 



Maximizing this overlap subject to a maximum feedback 
strength A yields the feedback coefficients 



A, = Asgn([nc,F,]), 



(43) 



This control solution, as for the bit-flip code, requires a 
controller to compute the feedback ( ^3|) . A natural ques- 
tion to ask is how the scaling of the classical computation 
behaves. In the appendix we show that the evolution of 
(2""*^)^ parameters must be calculated in order to com- 
pute the feedback for an [[n, k, d]] code, which at first 



does not seem promising. However, if one encodes mk 
qubits using m copies of an [[n, k, d]] code, as might well 
be the case for a quantum memory, the SME (^) will not 
couple the dynamics of the m logical qubits; and, as in 
the bit-flip case, the initial condition for the controller's 
integration can still be the completely mixed state in the 
total codespace. Then the relevant scaling for this sys- 
tem, the dependence on m, is linear: the number of pa- 
rameters is m(2"~'^')^. 



V. SIMULATION OF THE BIT-FLIP CODE 

In this section, we present the results of Monte Carlo 
simulations of the implementation of the protocol de- 
scribed in Section IV for the bit-flip code. 



A. Simulation details 

Because the bit-flip code feedback control scheme (|3l|~ 
^) uses a nonlinear feedback Hamiltonian, numerical 
simulation is the most tractable route for its study. To 
obtain Pc{t), we directly integrated these equations us- 
ing a simple Euler integrator and a Gaussian random 
number generator. We found stable convergent solutions 
when we used a dimensionless time step jdt on the order 
of 10~^ and averaged over 10'' quantum trajectories. As 
a benchmark, a typical run using these parameters took 
2-8 hours on a 400MHz Sun Uffra 2. We found that 
more sophisticated Milstein integrators converged 
more quickly but required too steep a reduction in time 
step to achieve the same level of stability. All of our sim- 
ulations began in the state Pc{0) = |0)(0| because it is 
maximally damaged by bit-flipping noise and therefore it 
yielded the most conservative results. 

We used two measures to assess the behavior of our bit- 
flip code feedback control scheme. The first measure we 
used is the codeword fidelity Fcw{t) = tr(pc(0)pc(i)), the 
overlap of the state with the target codeword. This mea- 
sure is appropriate when one cannot perform strong mea- 
surements and fast unitary operations, a realistic scenario 
for many physical systems. We compared Fcw{t) to the 
fidelities of one unprotected qubit Fi{t) = ^{1 + e~'^^*) 
and of three unprotected qubits -^^(t) = (Fi{t))^. 

The second measure we used is the correctable overlap 



where 



n(t) = tr {pc{t)Ilcorr), 



po + XIIpoXII 

-IXIpoIXI + IIXpoIIX 



(44) 



(45) 



is the projector onto the states that can be corrected 
back to the original codeword by discrete quantum er- 
ror correction applied (once) at time t. This measure is 
appropriate when one can perform strong measurements 
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and fast unitary operations, but only at discrete time in- 
tervals of length t. We compared Fcorr{t) to the fidelity 
F3 (t) obtained when, instead of using our protocol up to 
time t, no correction was performed until the final dis- 
crete quantum error correction at time t. As we showed 
in equation the expression for F^{t) may be calcu- 
lated analytically; it is F^it) = j{2 + Se-^T* - e'^''^) - 
l-372<2. 



B. Results 

We find that both our optimized estimate feedback 
scheme (|38) and our heuristically-motivated feedback 
scheme (^) effectively protect a qubit from bit-flip 
decoherence. In Figs. ^ and || we show how these 
schemes behave for the (scaled) measurement and feed- 
back strengths k/j = 64, A/7 = 128 when averaged over 
10^ quantum trajectories. Using our first measure, we see 
that at very short times, both schemes have codeword fi- 
delities Fcw{t) that follow the three-qubit fidelity F3{t) 
closely. For both schemes, Fcw{t) improves and surpasses 
the fidelity of a single unprotected qubit Fi{t). Indeed, 
perhaps the most exciting feature of these figures is that 
eventually Fcw{t) surpasses i^3(i), the fidelity achievable 
by discrete quantum error correction applied at time t. 
In other words, our scheme alone outperforms discrete 
quantum error correction alone if the time between cor- 
rections is sufficiently long. 




0.1 0.15 
Time 



0.25 



FIG. 2. Behavior of our protocol with optimized feedback 
(KSl) for parameters K/7 — 64, A/7 = 128, averaged over 
10 quantum trajectories. The analytical curves shown are 
as follows: the dashed line is the fidelity of one decohering 
qubit, Flit); the dashed-dotted line is the fidelity of three 
decohering qubits, Fz(t); and the dotted fine is the fidelity of 
an encoded qubit after one round of discrete error correction, 
F^j{t). Our simulation results are as follows: the solid line is 
the codeword fidefity Fcw{t), and the thick solid line is the 
correctable overlap Fcorr{t). 



Looking at our second measure in Figs. || and ||, we 
see that Fcorr{t) is as good as or surpasses F^{t) almost 
everywhere. For times even as short as a tenth of a de- 
coherence time, the effect of using our protocol between 
discrete quantum error correction cycles is quite notice- 
able. This improvement suggests that, even when one 
can approximate discrete quantum error correction but 
only apply it every so often, it pays to use our protocol 
in between corrections. Therefore, our protocol offers a 
means of improving the fidelity of a quantum memory 
even after the system has been isolated as well as possi- 
ble and discrete quantum error correction is applied as 
frequently as possible. 

There is a small time range from t = 0.01 to t ^ 0.05 
for the parameters used in Fig. || in which using our pro- 
tocol before discrete quantum error correction actually 
underperforms not doing anything before the correction. 
Our simulations suggest that the reason for this narrow 
window of deficiency is that, in the absence of our proto- 
col, it is possible to have two errors on a qubit (e.g., two 
bit flips) that cancel each other out before discrete quan- 
tum error correction is performed. In contrast, our pro- 
tocol will immediately start to correct for the first error 
before the second one happens, so we lose the advantage 
of this sort of cancellation. This view is supported by the 
fact that Fcorr (t) in our simulations always lies above the 
fidelity line obtained by subtracting such fortuitous can- 
cellations from F^{t). In any case, this window can be 
made arbitrarily small and pushed arbitrarily close to the 
beginning of our protocol by increasing the measurement 
strength k and the feedback strength A. 




0.1 0.15 
Time 



0.25 



FIG. 3. Behavior of our protocol with non-optimized feed- 
back ( p^ for parameters k/j = 64, A/7 = 128, averaged over 
10^ quantum trajectories. As in Fig. ^, the dashed line is 
Fi{t), the dashed-dotted fine is F^it), the dotted line is F^{t), 
the solid line is Fcwit) and the thick solid line is Fcw{t). Note 
that this feedback is qualitatively similar to that in Fig. ^ but 
does not perform as well. 

In Figs. ^ and ||, the Fcw{t) line is much more jagged 
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than the Fcorr{t) line. The jaggedness in both of these 
hnes is due to statistical noise in our simulation and is 
reduced when we average over more than 10"' trajecto- 
ries. The reason for the reduced noise in the Fcorrit) 
line has to do with the properties of discrete quantum 
error correction — on average, neighboring states get cor- 
rected back to the same state by discrete quantum error 
correction, so noise fluctuations become smoothed out. 

The improvement our optimized estimate feedback 
protocol yields beyond our heuristically-motivated feed- 
back protocol is more noticeable in Fcw{t) than in 

Fcorr 

(t) as seen in Figs. || and |[ Our optimized protocol 
acts to minimize the distance between the current state 
and the codespace, not between the current state and the 
space of states correctable back to the original codeword, 
so this observation is perhaps not surprising. In fact, op- 
timizing feedback relative to Fcorr{t) is not even possible 
without knowing the codeword being protected. Never- 
theless, our optimized protocol does perform better, so 
henceforth we shall restrict our to discussion to it. 

We investigated how our protocol behaved when the 
scaled measurement strength k/j and feedback strength 
A/7 were varied using the two measures described in Sec. 
Our first measure, the codeword fidelity Fcw{t), 



VA 



crosses the unprotected qubit fidelity Fi{t) at various 
times r as depicted in Fig. ^. This time is of interest be- 
cause it is the time after which our optimized protocol im- 
proves the fidelity of a qubit beyond what it would have 
been if it were left to itself. Increasing the scaled feed- 
back strength A/7 improves our scheme and reduces t, 
but the dependence on the scaled measurement strength 
k/j is not so obvious from Fig. ^. 




120 120 



k/ Y 



X/y 



FIG. 4. Time r at which Fcw{t) = ^1(1") as a function of 
measurement strength n/j and feedback strength A/7. This 
crossing time is the time after which our optimized protocol 
improves the fidelity of a qubit beyond what it would have 
been if it were left to itself. 

By looking at cross-sections of Fig. ^, such as at 



A/7 = 80 as in Fig. |5[ we see that for a given scaled feed- 
back strength A/7 there is a minimum crossing time r as 
a function of measurement strength k/j. In other words, 
there is an optimal choice of measurement strength K/7. 
This optimal choice arises because syndrome measure- 
ments, which localize states near error subspaces, com- 
pete with Hamiltonian correction operations, which co- 
herently rotate states between the nontrivial error sub- 
spaces to the trivial error subspace. This phenomenon 
is a feature of our continuous-time protocol that is not 
present in discrete quantum error correction; in the for- 
mer, measurement and correction are simultaneous, while 
in the latter, measurement and correction are separate 
non-interfering processes. 
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40 
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FIG. 5. Time r at which Fcw{t) = Fi{t) as a function of 
measurement strength k/7, keeping correction strength fixed 
at A/7 = 80. 

In order to study how our second measure, the cor- 
rectable overlap Fcorrit), varies with k and A, we found 
it instructive to examine its behavior at a particular time. 
In Fig. ||we plot Fcorr{t), evaluated at the time t = O.2/7, 
as a function of k and A. As we found with the crossing 
time T, increasing A always improves performance, but 
increasing k does not because measurement can compete 
with correction. Since ^3(6. 2/7) = 0.927, for all the k 
and A plotted in Fig. ^, using our protocol between dis- 
crete quantum error correction intervals of time O.2/7 
improves the reliability of the encoded data. 
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ment strength K/7 and feedback strength A/7. This quan- 
tity corresponds to the fidelity of a state given continuous 
error-correction up to 'yt — 0.2, at which point discrete er- 
ror-correction is performed. 



VI. CONCLUSION 

In many realistic quantum computing architectures, 
weak measurements and Hamiltonian operations are 
likely to be the tools available to protect quantum states 
from decoherence. Moreover, even quantum systems in 
which strong measurements and fast operations are well- 
approximated, such as ion traps ]34| , it is likely that 
these operations will only be possible at some maximum 
rate. Our protocol is able to continuously protect un- 
known quantum states using only weak measurements 
and Hamiltonian corrections and can improve the fidelity 
of quantum states beyond rate-limited quantum error 
correction. In addition, because our protocol responds 
to the entire measurement record and not to instanta- 
neous measurement results, it will not propagate errors 
badly and therefore has a limited inherent fault-tolerance 
that ordinary quantum error correction does not. 

We expect that our protocol will be applicable to 
other continuous-time quantum information processes, 
such as reliable state preparation and fault-tolerant quan- 
tum computation. We also expect that our approach 
will work when different continuous-time measurement 
tools are available, such as direct photodetection. Fi- 
nally, although current computing technology has limited 
our simulation investigation to few-qubit versions of our 
protocol, we are confident that many of the salient fea- 
tures we found in our thrce-qubit bit-flip code protocol 
will persist when our protocol is applied to larger codes. 
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APPENDIX A: FEEDBACK BASED ON THE 
COMPLETELY MIXED STATE 

Even though our quantum error correction feedback 
control scheme described in Section [V does not dis- 



tinguish between codewords, it is not obvious that we 
do not need to know the initial codeword to integrate 
its SME and calculate the relevant expectation values. 
Since we are interested in protecting unknown quan- 
tum states, this property is crucial to our scheme's suc- 
cess. Fortunately, for a large class of stabilizer codes, the 
computation of the feedback can be done by assuming 
the initial state is the completely mixed codespace state 
Pe — ^ Yl"J^i{I + Ml), which we prove here. 

We begin by defining the set G for the [[n, A:,(i]] code 
C with stabilizer S{C) as 

G = {as I a £ Pn,s G S{C), [s, a] = iff \a\ is even} , 

(Al) 

wher e \a\ denotes the weight of a as defined in section 

iiiBl. 



It is also useful to define the normalizer N(S) for the 
code as the group of operators which commute with every 
element in S{C). The elements of N{S)\S can be thought 
of as the encoded operations for the code — they move one 
codeword to another. 

We shall rewrite the conditions we require for the com- 
putation of the feedback to be insensitive to the initial 
codeword in terms of the Pauli basis coefficients Rg{p) 
which we define as follows. Let g = ® ... ® (Ji^, where 
ii . . .in take on the values x, y, z, I and cr/ = /. Then 



i?,(p)^tr(p.g)/2" = (.9)/2'' 



(A2) 



We can then formulate the problem in terms of proving 
conditions on G as follows: 

1. For every Rg used in our feedback scheme, g € G. 

2. For every g £ G and every pi and p2 in C, Rg{pi) = 
Rg{P2)- 

3. Evolution under the SME couples members of the 
set {Rg\g € G} only to each other. 

Theorem Let C be an [[n, 1,3]] ^ stabilizer code whose 
stabilizer S{C) has generators of only even weight and 



''The restriction to [[n, 1, 3]] codes is for simplicity of analy- 
sis; the proof may be extended to larger codes. Note that for 
an [[n, 1,3]] code, the Fi in the master equation ( ^ ) are all 
of the form o-j''\ where this notation denotes the weight-one 
Pauli operator ctj acting on qubit k. 
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whose encoded operations set N{S) \ S has elements of 
only odd weightlQ Then the conditions 1-3 above are sat- 
isfied; consequently, our scheme does not require knowl- 
edge of where the initial codeword lies in C. 
Proof: 

In this proof, any variable of the form is an ar- 
bitrary element of P„, and any variable of the form Sa 
is an arbitrary element of S{C). We prove each of the 
conditions listed above separately. 

Condition 1: By construction, G contains all M of 
the form M — Sia^l'\ where [si,aj'^^] ^ 0. These are 
precisely the operators used to compute the feedback in 
(^) for a code encoding one qubit. 

Condition 2: Let g = as G G and let p G C. We know 
either ae S,a e N{S)\S, or a ^ N{S). Suppose a& S. 
Then g G S acts trivially on all states in the codespace, 
so Rg — l/2"tr(pg) — 1/2" for this case. Now suppose 
a G ^{S) \ S. Then [a,s] ~ 0, and since as £ G, \a\ 
is even. But every element of N{S) \ S has odd weight 
by hypothesis, which is a contradiction. Hence a cannot 
be in N{S) \ S. Finally, suppose a ^ N{S). Then there 
exists some s' & S such that [a, s'] ^ 0; let s' be such an 
element. Then for \tp),\(l)) e C, 

(^|a|</.) = (V|as'|0> =-(V'|s'a|<^> 

= -(^|«|0)=O. (A3) 



weight of ak) 



N 



akSiSj 



OLkSr. 



so the con- 



dition above that [sj, iV] = becomes [sj,akSiSj\ = 
{ak[sj,sis.j\ + [sj,ak\siSj) [sj,ak] = 0. Therefore, 
[ak,Sm] = si[ak,Sj] -\- [ak,si]sj = [ak,si]sj which is zero 
or not depending on the original weight of ak- So if 
M — akSi is such that M G G, N — aks„i must fulfill 
that same condition, implying that N G G also. 

Similarly, Rn contributes to dRM through the master 
equation term [a^'^^ , p] when Na^''^ — M and [u^'"-' , N] ^ 
0. Now, M g G so M = aiSm, again with the appro- 
priate restriction on [a/ , Sm] depending on the weight of 
a/. Then N = o^j'^'^aiSm = ctnSm, so the condition above 
that {<yj''\ N} =/= becomes 

r (fe) (fe) 1 (fe)r (fe) 1 I r C^) i 



(fe) r (k) 1 IKl r IKI i 

= 0. (A4) 

We can now divide the analysis of this term into two 
cases. Case 1 occurs when af^ai has weight \ai\, im- 



plying that {ai,af^} 



0. Then {af\(jf\is,n} 



(fc) r (fc 
(fc) 



[c:- , s„i] — 0, which implies that [s,„,q;„] = 



(fc)r 
-a) >\s,, 



\Sm,(y] 'jaz+crj '[s,-n,ai\ = aj '[Sm,ai]. So [sm,a„] =0 



Hence for this case Rg 
these expressions for Rg must be the same no matter 
where p is in the codespace; therefore, for every g & G 
and pi,p2 e C, Rg{pi) = Rg{p2)- 

Condition 3: We prove this by considering dRj^j , where 
M e G: we will show that dRu = fi{RN\N € G}) 
for some real function /. Now, for any Ms Pn, 
dRM — Tr((i/5 M), where dp is given by the master equa- 
tion (42), and we can show condition 3 for each term of 
the master equation separately. First, substituting in the 
master equation shows that any term of the form 'D[c\pdt 
contributes either or the simple exponential damping 
term —2Rm to di?M if M and c commute or anticom- 
mute, respectively. 

As for the master equation term T-l[sj]dWjp, by writ- 
ing the master equation in the Pauli basis we can see 
that Rn contributes to dRjxj through this term pre- 
cisely when Nsj = M and {sj,N} ^ 0. Here we 
know that M g G, so we may write M = a^s; (with 
the appropriate restriction on [offcjS;] depending on the 



just when [s^jQ^/] — 0, which means that N G G since 
l/2"tr(pas) = 0. Note that |q;„| = |a/|. 

In Case 2, ct^'^^q;; has weight |q;; ± 1| 



[ai,cr'^'''>] 0. Then (|Aj becomes {cr)"' ,(7)"' ais^} 



(fe) ^(fe). 



,(fe). 



-3 

r (fe)i I (fe)r 

just when {s™, ai\ 
|a„| = la/ ±11. □ 



Sm\ = 0, which implies that [sm,a,i] = 

Oil} = (7f\sm,ai}. So [s„j,a„] = 
0, which means that N ^ G since 



^It is possible that this restriction may be able to be re- 
laxed; however, it is sufficiently general that it holds for the 
most well-known codes, including the bit-flip code, the five- 
bit code, the Steane code, and the nine-bit Shor code. This 
condition also ensures that G is consistent, i.e., if OjSk £ G 
and Qj — anSm, then a,i and SmSk also fulfill the conditions 
for an{smSk) to be in G. 



Thus we have shown the three conditions that all the 
R's used to compute the feedback are of the form Rni^g] 
that for a given M G G, Rm will be the same for any 
state in the codespace; and that evolution via the master 
equation mixes the i?'s of the form Rneg only with each 
other. Therefore, we can conclude that taking the initial 
state to be any state in the codespace, including the true 
initial state and the entirely mixed state, produces the 
same expression for the feedback when the master equa- 
tion is evolved conditioned on a measurement record, and 
so we do not have to know the true initial state to use 
our protocol. 

Another consequence of using the completely mixed 
state for feedback arises from the fact that doing so cor- 
responds to discarding information about the state of 
the system. Therefore, this procedure should reduce the 
number of parameters needed to compute the feedback. 
Unfortunately, this only leads to a modest reduction in 
the number of parameters, which can be found by using 
a simple counting argument. There are 2"/2''' = 2""'^ 
different error subspaces, including the no-error (code) 
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space, and if wc start with the completely mixed state 
in the codespace we do not need to worry at all about 
any movement within any of these spaces. Wo must only 
worry about which error space we are actually in, along 
with coherences between these spaces, so we find that 
^2"-fc)2 parameters are needed to describe the system. 
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